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ABSTRACT 


The  mixed-mode,  elasto-dynamlc  state  of  stresses  In  the  neighborhood  of  a 
rapidly  running  crack  tip  has  been  used  to  develop  a  relation  between  the 
Isochromatlc  fringe  orders  Ni)  and  Its  position  parameters.^  and  0.)  The  maximum 
shear  stress  Is  expressed  In  terms  of  the  stress  Intensity  factors  (1L,  «jj, 


and  other  higher  order  terms  Involving  the  mixed-mode  loading  for  a  crack 
propagating  at  constant  velocity.  A  graphics  package  based  on  these  derivation 
was  developed  for  mapping  the  Isochranatlcs  In  the  vicinity  of  a  running  crack 
tip  and  was  used  to  Illustrate  typical  mixed-mode  Isochromatics.  ^JThe^^ 
unsymmetry  associated  with  higher  order  terms  of  mixed-mode  stress  field  with 
the  mode  I  singular  stress  field  and  with/without  the  mode  II  singular  stress 


field  also  Is  Investigated.  Error  estimates  due  to  the  assumed  presence  of 

In  »  Kj  stress  field  was  found  to  be  significant  when  the  distance  from  the 
crack  tip  Is  more  than  4mm. 


*  University  of  Washington,  Department  of  Mechanical  Engineering,  Seattle, 
Washington  98195 
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INTRODUCTION 


While  much  has  been  discussed  theoretically  and  experimentally  on 
mixed-mode  crack  propagation  Cl— 43 *  all  discussions  Involved  only  quasi-static 
mixed-mode  crack  propagation.  As  for  rapidly  propagating  cracks*  there  exists 
ample  experimental  evidence  which  shows  that  dynamic  effects  significantly 
affect  the  crack  tip  stress  fields  [5-10].  Dynamic  analysis  of  a  rapidly 
propagating  crack  Is  also  required  for  Interpreting  the  observed  crack  curving 
and  crack  bifurcations  In  brittle  materials.  The  directional  stability  of  a 
rapidly  running  crack  has  been  shown  to  be  Influenced  by  crack  velocity  and  the 
non-singular  stress*  aQX,  which  Is  the  second  order  term  In  the  crack-tip 
stress  field.  This  stability  problem  was  discussed  theoretically  and  verified 
experimentally  with  dynamic  photoelastic  results  by  the  authors  [11*12]  and 
Rossmanlth  [13].  These  photoelastic  results  are  based  on  data  measurements  In 
the  vicinity  of  the  crack  tip  and  can  Introduce  errors  In  estimated  values  of 
fracture  parameters  due  to  crack  tip  blunting*  Insufficient  fringe  resolution 
and  Inaccurate  determination  of  crack  tip  location.  On  the  other  hand*  other 
far  field  errors  may  be  Introduced  when  data  away  from  this  crack  tip  zone 
together  with  higher  order  terms  In  the  crack  tip  stress  field  Is  used  to 
characterize  the  near  crack  tip  stress  field. 

Two  and  three  parameter  characterization  of  static  and  dynamic  stress 
fields  under  pure  mode  I  loading  condition  has  been  studied  exhaustively  by 
Rossmanlth  and  Irwin  [7].  The  effects  of  the  first  higher  order  term*  oqx, 
which  Is  the  lowest  order  non-singular  term  and  represents  the  stress 
biaxial Ity  at  the  crack  tip*  were  studied  by  Kobayashl  and  Ramulu  [9-12]  for 
static  and  dynamic  mixed-mode  conditions.  In  a  recent  Investigation* 
Rossmanlth  [14]  observed  that  the  unsymmetry  of  Isochromatics  can  be  generated 
by  superimposing  higher  order  terms  of  mode  II  stresses  on  mode  I  stress  field 


away  fraa  the  crack  tip.  Our  dynamic  photoelastic  analyses  of  crack  curving 
and  crack  branching  experiments*  however*  showed  that  the  unsymmetry  of 
Isochromatlcs  at  a  distance  as  close  as  1  to  3mm  from  the  crack  tip  was  always 
associated  with  dynamic  Kjj  and  hence  with  directionally  unstable  cracks. 
Because  the  relative  magnitude  of  such  dynamic  Kjj  small  with  respect  to  the 
dominating  dynamic  Kj  values  for  many  crack  propagation  problems  and  Its 
Influence  on  the  far  field  Isochromatlcs  Is  usually  small.  On  the  other  hand* 
unsymmetry  In  the  far  field  Isochromatlcs  may  not  require  the  presence  of 
dynamic  Kjj,  as  reported  by  Rossmanlth  [143  and  thus*  error  estimates  In  the 
fracture  parameters  due  to  the  assumed  presence  of  KII  In  a  Kj  stress  field  Is 
needed.  A  numerical  experimentation  based  on  the  theoretical  dynamic 
Isochromatlc  field  could  provide  not  only  the  errors  Involved  In  wrongful 
presence  of  dynamic  but  also  could  provide  Information  the  minimum  number 
of  higher  order  terms  needed  In  the  crack  tip  stress  field  and  region  of 
Influence  of  the  near  crack  tip  field  under  dynamic  conditions  for  accurate 
assessment  of  the  dynamic  crack  tip  stress  field. 

The  objective  of  the  present  Investigation  Is  to  study  the  effects  of 
higher  order  terms  on  the  mixed-mode  crack  tip  stress  pattern  of  a  rapidly 
running  crack  with  an  attempt  to  rationally  delineate  the  optimum  region  of 
data  collection.  The  effects  of  higher  crack  velocity  on  the  crack  tip  stress 
pattern  which  has  not  been  reported  previously  under  mixed-mode  loading  also 
are  presented. 

THEORY 

The  general  near  crack  tip  dynamic  state  of  stresses  under  mixed-mode 
conditions  was  given  recently  by  Nlshloka  and  Atlurl  [15]  In  terms  of  the  local 
rectangular  (x*y)  and  polar  (r*8)  coordinates.  The  three  rectangular  stress 
components  under  mode  I  and  mode  II  conditions  are  given  as 
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r  and  0  are- the  polar  coordinates  with  origin  at  the  moving  crack  tip  and  C, 

Cv  and  C2  are  the  crack  velocity,  dilatations!  and  dlstortlonal  wave 
velocities,  respectively. 

The  general  solution  expressed  in  Equation  Cl)  yields  the  singular  stresses 

when  n  =  lj  l.e.,  Aj^  =  Kj  and  Ajj^  =  Kjj,  which  are  stress  intensity  factors 
of  mode  I  and  mode  II,  respectively.  The  constant  stress,  °ox*  is  related  to 
the  higher  order  terms  (HOT)  for  n  >  2  as 
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The  n  =  2  term  exists  only  In  stress  component  of  mode  I  loading  and  is 

zero  In  all  other  exponents  of  mode  I  and  .  ode  JI  stress  fields..  The 
Influence  of  these  and  third  order  terms,  n  =  3,  on  the  shape  of  the 
Isochromatics  surrounding  the  crack  tip,  are  studied. 

The  analyses  of  mixed-mode  Isochromatics  Is  made  by  computing  the  maximum 
shear  stress  contours  around  the  crack  tip.  The  maximum,  in-plane  shear 
stress,  t^,  is  related  to  the  cartesian  stress  components  as 


<2r  )2  =  Co  -o  )2  +  2t  2 
V  yy  xx'  xy 


(3) 


The  analytical  expression  for  mixed-mode  isochromatics  in  the  vicinity  of  the 
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crack  tip  is  obtained  by  combining  the  maximum  shear  stress*  r  ,  and  the  stress 

m 

optic  law. 
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where  N  is  the  fringe „ order#  fa  is  the  stress-optical  constant*  and  h  is  the 
specimen  thickness.  The  combination  of  Equations  (1-4)  yields  a  highly  complex 
relation  which  can  be  written  In  a  functional  form  as 

N  =  F(f(j*h t  r#  Kj,  Kj j#aQX,  C»  C^»  C2*  ®#  HOT)  C5) 

The  Isochranatlc  field  near  the  crack  tip  for  given  fracture  parameters  of 

Kp  K11#  °ox»  and  the  model  parameters  of  foand  h  as  well  as  the  crack 
velocity*  C»  may  be  generated  by  computing  N  at  a  large  number  of  special 
points.  Equation  (5)  Is  programmed  to  plot  the  dynamic  lsochromatlcs  where  the 
computer  routine  Involved  Is  similar  to  that  discussed  In  References  [33  and 
C93 .  -  The  values  were  computed  for  an  array  of  points  within  a  region  of  ± 
25mm  with  the  crack  tip  In  the  center. 


KFAR  CRACK  TIP  ISQCHRQMAT.ICS 

As  an  accuracy  check  of  Equation  (5)»  the  dynamic  lsochromatlcs  should 
coincide  with  that  of  References  [3]  and  [93  for  the  limiting  case  of  crack 
velocity  C  ->  0  as  well  as  C  =  0.15Cp  respectively*  for  Identical  values  of 
Kp  Kjp  and  oQX.  These  fracture  parameters  Kj*  Kjp  and  the  other  higher 
order  terras  through  the  third  order  were  varied  and  crack  tip  Isocbromatlc 
fringes  were  constructed.  The  size*  shape*  and  orientation  of  these  fringe 
patterns  depend  sirongly  on  the  combination  of  parameters  chosen  [3*4*93.  In 
order  to  visualize  the  effects  of  the  higher  order  terms*  theoretical 
lsochromatlcs  were  generated  using  the  same  model -fringe  constant  of  fc/h  » 
1.73  MPa/fringe  <250psi/fr1nge)*  Kj  0.876  MPa  fm  (800psf  ^[7n>  and  Kn  of 
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Reference  [33.  A  constant  crack  velocity  of  C  -  0.20C^  and  0.30C^  was  chosen 
to  generate  the  dynamic  Isochromatics.  In  addition*  «  2, 400m/ s  (94*300 
In/s)  and  C2  «  1,160  m/s  (45*800  1n/s)»  Kjj/Kj  >  0  were  assumed  In  the  dynamic 
analysis.  Static  Isochromatics  were  generated  by  setting  C  «  0.01  in 
Equation  (5)  and  these  plots  were  In  agreement  with  that  of  References  [3]  and 
[4] . 

The  static  and  dynamic  Isochromatics  listed  In  Figures  1  through  8  can  be 
broadly  classified  as: 

a)  Mode  I  (Ajn  >  o  *  0#  Ajj  n*  0) 

b)  Mode  II  (AIn  .  o,  AIIn  ^  0) 

c)  Mixed  Mode  (^IIl/Ajj  »  kII/Kj  ■  0.25  and  AIn  ■  Ajjn  >  0  when  n  >.  2) 

The  higher  order  terms  A^  ^  0  accounts  for  the  stress  blaxlallty  under 

mode  I  loading  condition  and  Is  Independent  of  the  radial  distance  and  the 

angle.  The  Aj3  and  Ajjj  terms  represents  the  Influence  of  the  near  boundaries 
on  the  Isochromatics  and  Its  magnitude  and  sign  depends  on  the  loading 
configuration  [43.  In  the  following  discussion*  seme  representative  static  and 
dynamic  Isochromatics  In  the  vicinity  of  the  crack  tip  are  shown  In  Figure  1-8 
for  mode  I*  II*  and  mixed-mode  loadings. 

PURE  MODE  I  (Ajn^  0;  Ajjn  ■  0  for  n  i  2) 

The  analytically  generated  pure  mode  I  crack  tip  stress  patterns  with  the 
effect  of  higher  order  terms  of  n  £  2*3*  which  corresponds  to  r®  and  r1^* 
respectively,  for  a  crack  velocity  of  C/Cj  «  o.2*  Is  shown  In  Figure  1.  The 
fringe  pattern  tilt  forward  for  AJ2  <o  and  backwards  for  AI2  >0.  The  added 
third  order  term  of  Aj3  <o  increases  the  maximum  apogee  distance  rmax  of  the 
Isochromatlc  for  Aj2  <0  and  decreases  r,,^  when  Aj2  >0.  The  effect  of  A^  term 
Is  noted  as  the  distance  from  the  crack  tip  Increases  and  appears  to  be  an 


7 


essential  term  In  the  evaluation  of  fracture  parameters  using  far-fleld  data. 
Detailed  discussion  of  mode  i  Isochromatics  can  be  found  In  Reference  [7]  and 
[9] . 

PURE  MODE  II  (AIn  «  0;  AIIn  £  0  for  n  >  2) 

Figure  2  shows  the  pure  mode  II  static  Isochromatlc  stress  field  for 

AII1  f  0  and  ^113  >®*  The  effect  of  Ajj3  term  Is  significant.  When  Aju  >  0 
and  AII3  *  0»  the  Isochromatics  are  symmetric  with  respect  to  x  and  y.  When 
AII3  1s  added  to  the  pure  mode  II  singular  stress  field*  only  the  x  axis 
symmetry  Is  maintained. 

Figure  3  shows  the  pure  mode  II  singular  crack  tip  stress  patterns  (n  *  11 
with  Increasing  crack  velocity  C/C^#  As  the  crack  velocity  Increases*  the  size 
of  the  Isochromatlcs  ahead  of  the  crack  tip  decreases.  Although  there  Is  no 
experimental  evidence  of  a  crack  propagating  rapidly  under  pure  mode  II 
loading*  this  numerical  experiment  demonstrates  that  the  crack  velocity 
significantly  affects  the  Isochromatlc  patterns. 

MIXED  MODE  (AJn  >  0,  AIIn  t  0) 

Superposition  of  a  slight  shearing  stress  component  with  and  without  the 
singular  stress  component  destroys  the  symmetry  of  the  mode  I  Isochromatlc 
pattern  [14].  The  mixed-mode  cases  are  classified  Into  three  groups: 

a)  Mixed-mode  with  dominant  mode  I  and  higher  order  terms. 

b)  Mixed-mode  with  dominant  mode  II  and  higher  order  terms. 

c)  Unsymmetrlc  Isochromatlcs  associated  with  higher  orders  of  mode  II 
terms  In  (a)  and  with  higher  orders  of  mode  I  terms  In  (b). 

In  the  following*  some  representative  static  and  dynamic  Isochromatlcs  are 
presented. 


Figure  4  shows  the  pure  mixed-mode  static  and  dynamic  singular  stress 

patterns  for  kjj/Kj  »  0.25.  The  fringe  patterns  shown  are  symmetric  and  the 
axis  of  symmetry  rotated  clockwise  about  25  degrees  for  *  o.25»  as 

observed  by  Rosamanlth  [4].  For  Kjj/Kj  >  0.25#  the  axis  of  symmetry  rotates 
counterclockwise  about  25  degrees.  Furthermore*  the  isochromatlc  fringe 
pattern  distorts  at  higher  crack  velocities  by  shortening  and  stretching  of  the 
upper  loops  and  lower  loops#  respectively.  For  C/Cj  .  o.30#  the  upper  fringe 
loops  tilt  counterclockwise  while  the  lower  fringe  loops  do  not  rotate  any 
further. 

Figure  5  shows  the  mixed-mode  crack  tip  stress  pattern  with  a  dominant  mode 
I  stress  field.  The  Isochromatics  Increase  In  size  for  the  positive  higher 
order  terms  of  AI2,  aI3#  and  for  crack  velocities  of  0.2^  and  0.3^.  The 
upper  loop  again  tilts  counterclockwise  and  flattens  as  the  magnitude  of  C/Cj 
Increases  to  0.3#  as  shown  In  Figure  5a.  Similar  effects  were  noted  for  the 
experimentally  observed  crack  velocities  of  C/Cj^  *  o.2»  with  the  addition  of 
the  negative  higher  order  terms  of  AJ2  and  AI3.  Addition  of  the  third  higher 
order  term#  Ajj,  regardless  of  the  sign  of  the  A^  In  the  stress  field# 
Increased  the  fringe  size.  Velocity  effect  Is  significant  at  C/C ^  *  o.30  for 
added  terms  of  Aj2  >  o. 

The  Influence  of  Ajj  and  Ajq  terms*  which  represent  the  boundary  Influence 
on  the  Isochromatics#  Is  shown  In  Figure  5b  under  mixed-mode  loading#  with  AI2 
*  0.  For  the  equal  magnitudes  of  Ajj  and  Ajj3#  the  fringe  loop  size  Increased 
significantly  and  exhibited  severe  unsymmetrlc  stress  patterns  away  from  the 
crack  tip.  However*  In  the  near  vicinity  of  the  crack  tip#  the  fringe  pattern 
shape  and  orientation  Is  quite  similar  to  the  fringe  patterns  of  Figure  4. 


This  demonstrates  the  fact  higher  order  terms  must  be  Included  In  the  stress 
field  characterization  when  far-fleld  fracture  data  are  used. 

Mode  II  Dominant  Loading 

The  crack  velocity  not  only  affects  the  mixed-mode  Isochromatics  for  mode  I 
dominant  loading  but  also  for  mode  II  dominant  loading.  Figure  6  shows  the 
theoretically  generated  Isochromatics  with  Increasing  mode  II  loading 
representing  mode  II  dominant  Isochromatlcs.  Under  loading*  the  fringe 
loops  do  not  terminate  at  the  crack  tip  and  the  fringe  patterns  become 
Increasingly  dissimilar.  Increases  In  mode  II  loading  rotates  the  fringe 
pattern  counterclockwise  with  changes  In  shape  and  size. 

Details  of  the  changes  In  the  shapes  of  the  static  mixed-mode  Isochromatlcs 
are  discussed  In  Reference  [4]  and*  thus*  will  not  be  reproduced  here.  In 
essence*  the  corresponding  dynamic  Isochromatlcs  are  found  to  be  larger  and 
essentially  follow  the  general  shape  of  the  static  Isochromatlcs.  Similar 
results  are  shown  In  Figures  1  through  9  for  the  experimentally  observed  crack 
C/Cx  *  0.2. 

Usymmetrlc  Isochromatlcs 

Usymmetrlc  Isochromatlc  patterns  can  be  generated  by  adding  AJI3  >  o*  with 
and  without  Aj3  >  o.  However*  significant  change  In  Isochromatlc  patterns  are 
noted  when  the  sign  of  the  remote-stress  component  oqx,  i.e. »  Aj2»  Is  changed. 
As  shown  In  Figure  5*  the  fringe  loops  In  the  angular  range  of  0  to  180° 
decrease  In  size  for  a  positive  aQx  and  mode  II  stress  Intensity  factor.  The 
higher  order  terms  (n*3)  appear  to  have  a  negligible  effect  on  the 
Isochromatlcs  In  the  vicinity  of  the  crack  tip.  The  loss  of  symmetry  and  the 
distortion  of  the  Isochromatlcs*  however*  become  more  pronounced  when  higher 
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order  terms  without  singular  term  mixed-mode  loading  (Ajn  4  0*  Kjj  =  0,  AIIn  4 
0)  are  added.  The  fringe  loops  ahead  of  crack  tip  also  do  not  terminate  at  the 
crack  tip. 

Figure  7  shows  the  effect  of  the  higher  order  terms  of  mode  II 
Isochramatlcs  with  the  mode  I  stress  field.  When  Ajjj  >  0  Is  superposed  onto 
the  mode  I  stress  field  with  Kj  >  g  an(j  Kjj  b  0,  -the  resulting  Isochromatlc 
stress  pattern  loses  Its  x-  and  y-axls  symmetries  for  all  Aj2,  Increase  In 
crack  velocity  distorts  this  Isochromatlc  pattern.  Similar  changes  also  are 
seen  In  the  case  of  Figure  8.  when  the  higher  order  terms  of  mode  I  stress 
field  with  Kj  s  0,  Aj2  >0.  and  Ajj  »  Ajjj  £  0  are  superposed  onto  pure  mode  II 
singular  stress  field. 

Slight  unsymmetry  In  load  are  often  observed  experimentally  and  are  modeled 

with  Kj j  *  g,  Ajj3  <  0  In  the  presence  of  aQx.  Such  loading  Is  superposed  onto 
the  mode  I  crack  shown  In  Figure  9a.  A  genuine  mixed-mode  crack  tip 

Isochromatics  Is  modeled  with  Kjj/Kj  -  0.10  and  Is  also  shown  In  Figure  9b. 
Both  Isochromatics  exhibit  slight  unsymmetry  and  their  shapes,  sizes,  and 
orientation  are  Influenced  not  only  by  the  crack  velocity.  C/Cj,  but  also  by 
the  ratio  of  mixed-mode  stress  Intensity  factors. 

DISCUSSIONS 

Figures  1.  3.  5.  and  6  show  that  the  crack  velocity  affects  the 
Isochromatics  associated  with  mode  I.  II.  and  mixed-mode  loading,  and  also  In 
mode  II  dominant  loading.  Increasing  crack  speed  tends  to  reduce  the 
directional  stability,  which  Is  generated  by  compressive  stress  parallel  to  the 
crack.  Irrespective  of  the  modes  of  deformation  as  noted  In  the  recently 
proposed  crack  curving  criteria  [11].  The  added  higher  order  terms  change  the 
shape  and  size,  which  are  governed  by  the  sign  of  these  terms,  of  the  transient 
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Isochromatics.  Although  mode  II  stress  Intensity  factor  Is  zero*  1.e.»  = 

0»  In  Figure  9a,  the  presence  of  other  higher  order  terms  of  the  opening  and 
shearing  mode  stress  field  generated  unsymmetrlc  Isochromatics.  Figure  9b 
shows  only  a  slight  unsymmetry  observed  In  the  presence  of  a  smaller  amount  of 

KII^KI*  $uch  unsymmetrlc  Isochromatics  are  typically  associated  with  curved 
crack  and  post  branching  cracks. 

In  order  to  demonstrate  the  errors  Involved  by  assuming  a  mixed-mode  crack 
tip  stress  field  <Kp  Kjp  and  °ox>»  the  grossly  unsymmetrlc  Isochromatics  of 
Figure  9a  Is  chosen  for  a  numerical  experiment.  This  numerical  experiment  will 
provide  assessments  of  the  effective  region  of  the  crack-tip  stress  field  and 
of  the  terms  needed  for  obtaining  accurate  fracture  parameters.  By  applying 
the  overdetermlnlstlc  procedure  [9,163,  the  fracture  parameters  with  added 
higher  order  terms  were  evaluated  using  the  fringe  orders  3,  4,  and  5  In 
Figure  9a.  A  total  of  20  data  points  which  were  located  In  the  vicinity  of  the 
maximum  radial  distance  from  the  crack  tip  In  each  fringe  order  were  used. 
Data  points  were  taken  over  the  range  2mm  <  rmax  <  8.5mm,  at  several  values  of 
r.  Table  1  shows  the  two-pa rameter  (KJf  ow),  three- pa raneter  (Kj,  Kn,  oox), 
and  five-parameter  (Ajn,  Ajjn  for  n  *  1,2,3),  which  were  determined  by  this 
numerical  experiment  and  the  exact  values  associated  with  the  theoretical 
Isochromatics.  These  results  show  that  the  fracture  parameters  generated  by 
this  far-fleld  data  result  In  an  apparent  KII  In  addition  to  Kj.  This  finding 
Is  In  agreement  with  the  results  of  Rossmanlth  [143,  where  Kjj  diminishes  as 
the  data  points  closer  to  the  crack  tip  are  used  In  data  reduction.  The 
estimated  fracture  parameters,  which  were  obtained  by  using  five-parameter, 
also  overestimated  the  stress  Intensity  factor.  Again,  this  error  decreased 
when  the  nearest  fringe  was  used.  These  unsymmetrlc  patterns  In  the  absence  of 

are  found  to  be  present  only  when  r  Is  greater  than  about  4mm. 
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w 

In  the  presence  of  a  small  mode  II  stress  Intensity  factor*  the 
analytically  generated  Isochranatlcs  exhibits  unsymmetry  at  a  radial  distance 
as  small  as  r  *  1mm  and  Is  always  associated  with  a  kinked  crack.  Table  1  also 
shows  that  the  five  parameter  method  appears  less  desirable  with  added 
difficulty  In  numerical  convergency.  Also*  the  constant  crack  speed  assumption 
Is  reasonable  only  If  the  measurement  region  Is  restricted  to  a  moderate  size 
as  specified  In  Reference  [173.  Finally*  this  numerical  experimentation 
confirms  the  general  validity  of  the  th ree-pa rameter  method. 

Typical  mixed-mode  experimental  and  theoretical  Isochranatlcs  are  shown  In 
Figures  10  through  12.  The  data  points  for  generating  the  theoretical 
Isochranatlcs  were  taken  over  a  region  of  radius  1  to  4mm*  centered  at  the 
crack  tip.  The  theoretically  generated  Isochranatlcs*  with  and  without  higher 
orders*  matched  the  experimental  fringe  patterns  within  the  sampled  region 
around  the  crack  tip. 

Figure  10  shows  one  enlarged  frame  out  of  a  16-frame  dynamic  photoelastic 
record  of  a  curved  crack  In  a  9.5mm  thick*  254  x  254mm  Homallte-100 
single-edge- notch  (SEN)  specimen  loaded  under  flxed-grlpped  tension  [11].  Kj, 

KH»  and  oox  were  estimated  by  the  three- pa  rameter  method  and  Is  used  to 
reconstruct  the  experimental  Isochromatics.  The  mixed-mode  effect  Is  visible 
for  a  radial  distance  as  small  as  1mm  and  the  theoretical  and  experimental 
Isochranatlcs  agree  reasonably  well  within  the  region  of  about  4mm.  The 
Isochranatlcs  associated  with  the  continuously  curving  post-branched  crack  In  a 
wedge-loaded  rectangular  double  cantelever  beam  (WL-RDCB)  specimen  [12]  Is 
shown  In  Figure  11.  The  noticeable  unsymmetry  Is  not  due  only  to  the  higher 
order  terms  of  mode  II  loading*  but  to  the  genuinely  mixed-mode  crack  tip 

deformation  of  Kjj/Kj  *  -0.17.  The  analytically  generated  Isochromatics 
matched  well  with  the  experimental  Isochranatlcs  within  the  sampled  region. 


-»v. 
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Figure  12  shows  the  arrested  branch  crack  with  a  dominant  mode  II  crack  tip 
deformation  during  unloading  process  In  a  3.2mm  thick*  127  x  225mm 
polycarbonate  single-edge  notch  specimen  [183.  The  theoretical  and 
experimental  Isochromatics  matched  better  within  a  crack  tip  zone  of  4mm  when 
the  three  parameter  method#  without  the  higher  order  terms#  was  used  for  data 
reduction.  Thus*  fracture  parameters  associated  with  directionally  stable  and 
unstable  cracks  can  be  determined  with  reasonable  accuracy  by  the  three 
parameter  method  as  long  as  the  data  points  are  within  a  circular  region  of 
radius  4mm  from  the  crack  tip. 

Finally#  the  possible  errors  In  crack  curving  angles  which  were  generated 
by  using  the  mixed-mode  three  parameter  (Kj,  Kjj#  °ox)  characterization  of 
unsymmetrlc  Isochromatics  of  Kj  field#  were  estimated  by  using  the  fracture 
parameters  generated  by  the  nearest  fringe  order  data  In  Table  1.  The 
estimated  fracture  parameters  with  the  nearest  or  highest  fringe  order  from 
Table  1  yield  an  apparent  mixed-mode  stress  Intensity  factor  ratio  of  Kjj/Kj  * 

-0.026.  By  assuming  a  critical  radius  of  rc  =  i.3mm  [11,123#  the  parameter 

°ox 

i~c— —  #  which  controls  the  direction  of  crack  propagation*  becomes  -0.0569. 
This  data  Is  used  to  predict  the  crack  curving  angles  by  using  the  maximum 
circumferential  stress  theory. 

Figure  13  shows  the  Influence  of  kjj/Kj  on  the  analytically  predicted  crack 

^  ^OX 

kinking  angles  for  ^rc  -jr-  =  -0.05690#  0.0#  and  0.0569.  A  mixed-mode  stress 
Intensity  factor  ratio  of  Kjj/Kj  =  -0.0155  or  -0.036  with  a  positive  or 

P  Oqv 

negative  4  c »  +  0.0569#  respectively#  Is  required  to  generate  a  measurable 
Ki  o 

crack  curving  angle  of  3  degrees.  Thus#  In  the  presence  of  a  negative^  - K°*» 

_  0OX 

the  magnitude  of  Kjj/Kj  must  double  that  required  for  a  posu„.-Fc  tf  <" 
order  to  generate  a  visible  crack  curving.  On  the  other  hand#  a  hypothetical 

Kjj  ■  -0.026  Kj  with  ^  °ox  =  -0.0569  will  result  In  a  crack  kinking  angle  of 

Ki 
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1.8  degrees*  which  Is  negligible.  By  chang  1  ngJr  ~^m  0.0569*  the  kinking 

~  C  Kt  „ 


angle  Increases  to  4.2  degrees*  which  Is  noticeable.  Such  positive’!?  — -  , 

c  Kj 

however*  Is  only  remotely  possible.  The  predicted  crack  curving  angle  from 
the  evaluated  fracture  data  of  the  unsymmetrlc  Isochromatics  Is  at  the  most  1.8 
degrees  and  Is  barely  noticeable.  Possible  errors  which  are  generated  by  using 
three  parameter  method  (Kj,  Kjj*  °ox)  In  evaluating  the  unsymmetrlc 
Isochranatlcs  without  a  singular  mode  II  Is  negligible  If  the  photoelastic  data 
Is  confined  to  a  circular  region  of  radius  4mm  from  the  crack  tip. 

The  above  results  address  the  question  of  the  needed  terms  In  the  dynamic 
mixed-mode  crack  tip  stress  field*  and  also  Indicate  the  limitation  In  crack 
tip  zone  size  for  accurate  estimation  of  stress  Intensity  factors.  In  order  to 
define  the  conditions  under  which  the  higher  order  terms  need  to  be  Included  In 
the  dynamic  stress  field  remains  a  topic  worthy  of  detailed  Investigation. 


CONCLUSIONS 

1.  The  higher  order  terms  significantly  affect  the  size  and  shapes  of 
dynamic  crack-tip  Isochromatics  generated  by  pure  mode  I*  mode  II*  and  under 
mixed-mode  loadings.  In  particular*  the  first  higher  order  term  or  remote 
stress  component  changes  the  crack  tip  stress  pattern  shape  and  size  In  all 
modes  of  loading. 

2.  The  higher  order  terms  of  n  -  3,  l.e.,  AJ3,  aII3*  significantly  affect 
the  Isochromatics  and  must  be  used  In  analyzing  the  pure  mode  II  fracture. 

3.  Crack  velocity  of  C/C^  >  o.2*  alter  significantly  the  size  and  shape  of 
the  crack  tip  stress  pattern  In  all  modes  of  loading. 

4.  The  higher  order  terms  (n  >  3)  have  little  Influence  on  the 
elastodynamlc  mixed-mode  Isochrcmatlcs  for  a  radial  distance  less  than  4mm. 

5.  The  Isochromatics  can  become  unsymmetrlc  In  the  presence  of  higher  order 
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terns  In  mode  II  stress  field  without  the  presence  of  the  singular  term 
Involving  KII- 

6.  By  using  the  three  parameter  model  one  can  estimate  reasonably  accurate 
fracture  parameters  associated  with  directionally  stable  and  unstable  cracks 
within  the  maximum  radial  distance  of  1  to  4mm. 
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TABLE  1 


Evaluation  of  Fracture  Parameters 
Associated  with  the  Fringe  Pattens  Shown  In  Figure  12a 

™  ~  ESTIMATION 

EXACT  1X0  PARAMETER  THREE  PARAMETER  FIVE  PARAMETER 

VALUES 

FRINGE  ORDER  345345345 

0.876  0.810  0.820  0.820  1.049  0.954  0.873  0.986  1.041  0.892 

(MPaJrn) 

AI2  -0.58  -1.20  -0.78  -0.72  -0.73  -0.50  -0.59  -0.83  -0.37  -0.62 
(MPa) 

AI3  0  . -  —  —  —  -0.69  -0.17  -0.35 

(MPaAflS) 

Kn  0  —  —  —  -0.043  -0.023  -0.019  0.077  -0.043  -0.022 

(MPatfn) 

Am  -0.28  —  —  —  —  —  --  0.35  0.30  -0.35 

<MPa/T«) 
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Figure  8.  Unsymmetr ic  laochromatics .  Mode  II  ( K j. „  =  0.876MPa4m)  Dominant  Loading,  Mode  I  Higher  Order 
Terms  and  Crack  Velocity  C/C  =  0.2 
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Experimentally  and  Theoretically  Generated  Isochromatics  Associated  With  a  Curved 
Crack  in  Single-Edged-Notch(SEN)  Specimen  Under  Fixed  Grip  Loading 
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Figure  13.  Crack  Kinking  Angles  Associated  With  Mode  I  Dominant  Loading 
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